Annalen der Physik, 21 March 2013 



Organic glasses: cluster structure of the random energy land- 
scape 

Sergey V. Novikov* 

A.N. Framkin Institute of Physical Chemistry and Electrochemistry, Leninsky prosp. 31, 119991 Moscow, 
Russia 

Received 1 September 2009, accepted 5 September 2009 

Key words Random medium, energetic disorder, clusters. 
PACS 05. 10. -a, 05.40.-a 

An appropriate model for the random energy landscape in organic glasses is a spatially correlated Gaussian 
field. We calculated the distribution of the average value of a Gaussian random field in a finite domain. The 
results of the calculation demonstrate a strong dependence of the width of the distribution on the spatial 
correlations of the field. Comparison with the simulation results for the distribution of the size of the 
cluster indicates that the distribution of an average field could serve as a useful tool for the estimation of 
the asymptotic behavior of the distribution of the size of the clusters for "deep" clusters where value of the 
field on each site is much greater than the rms disorder. We also demonstrate significant modification of the 
properties of energetic disorder in organic glasses at the vicinity of the electrode. 
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1 Introduction 

Most disordered organic materials are molecular glasses having very low concentration of intrinsic free 
carriers. They demonstrate measurable conductivity only after injection of carriers by the action of the 
laser pulse or strong electric field. At the same time, organic glasses usually have high concentration of 
dipoles and quadrupoles. Due to almost zero concentration of free carriers and lack of screening such 
molecules provide long range contribution to the overall energetic disorder for charge carriers. Long range 
sources inevitably lead to the strong spatial correlation of the random energy landscape U (r): in organic 
polar materials (dipolar glasses) correlation function of random energy C(r) = (U (r)U(Q)) decays as 1/r 
and in nonpolar materials (quadrupolar glasses) it decays as 1/r 3 JTJE1E1- In this paper the most attention 
will be paid to cluster structure of the dipolar glass. Using an exact analytic calculation as well as computer 
simulations it was shown that in the dipolar glass the random energy U(r) is a Gaussian random field if the 
average distance between dipoles is not significantly greater than the minimal distance a (a lattice scale) 

ma. 

Long range correlations mean that sites with close values of random energy tend to group together and 
form clusters (cluster is defined as a set of connected sites, where all of them have site energy U greater 
than some boundary energy Uq). Until recently very little has been known about cluster distribution on size 
n s in correlated Gaussian medium (here s is the number of sites in a cluster). We suggested an approximate 
method for the analytic estimation of n s for "deep" clusters with Uq er (here a is rms disorder and, by 
definition, C(0) = a 2 ). It turns out that a very simple calculation of the distribution of the average value 
of the correlated Gaussian field in finite domain containing s sites provides a very accurate estimation of 
the leading asymptotic for the distribution of large clusters with s 3> 1 and Uq 3> o. 

We also discuss briefly how the correlated properties of the energetic disorder in dipolar glasses are 
modified in the vicinity of a conducting electrode. 
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2 Distribution of domain averages for Gaussian energy landscapes 

The simplest domain average is the distribution Pv(Uq) of the average value Uq of the random energy 
U (f) in a domain with volume V (here we consider a spatial average and use the same notation Uq for the 
average energy) 

Pv(U ) = ( 6 (yf drU{f)f v (r) - U^j ^ , (1) 

where /y(r) equals 1 inside the domain and outside. For the Gaussian random field this distribution can 
be calculated in the exact way j<5| 

P V (U ) = -/== cxp (~^-) . K = J d? d?i fv (r)C(r - ri)/y (fi ) . (2) 

For the uncorrected Gaussian field C(r) oc S(r), hence, K oc V and does not depend on the shape of the 
domain. In the opposite case of the very strong spatial correlation, where C(r) m const at the microscopic 
scale, K oc V 2 and again does not depend on the shape of the domain. This extreme example looks quite 
unusual but, nonetheless, this is exactly the case for the Gaussian random energy landscape provided by 
the random uncorrected distribution of static charges if concentrations of charges of the opposite sings are 
equal 17] . In that case for the transport layer of the finite thickness L the only relevant spatial scale for 
the correlation function is L and K oc V 2 for V <C L 3 . We can conclude that the whole variety of the 
dependences of K(V) is bounded by these two limiting cases. Of course, in the typical situation K does 
depend on the shape of the domain, and these extreme cases are exceptional. 

The most studied case (apart from the trivial case of the uncorrelated Gaussian field) is the dipolar glass 
model, where C(r) w Aa 2 a/r with A as 0.76 (TJ. In that case for the spherical domain with radius Rq [6 | 

K = ^Aaa 2 R 5 oc V 5 / 3 . (3) 
1 5 

If a domain has an arbitrary shape but still could be characterized by a single linear scale I oc V 1 / 3 , then 
K oc V 5 / 3 just because of dimensionality argument, though the coefficient of proportionality depends on 
the actual shape of the domain. 

Importance of the distribution Py {Uq) stems from the fact that it can serve as a good estimation for the 
number of true clusters if Uq 3> o. Indeed, the number of spherical domains ny{Uo) per unit volume, 
having an average energy greater than Uq, is approximately equal to 

u v (Uq) xs v j ua dUP v {U) = -erfc (^=j ► JJ -^ J = cx P , (4) 

here the coefficient 1/V reflects the number of non-overlapping independent domains in a finite sample. 
One may expect that (O gives a reasonable estimation for the number n s of the true clusters, i.e. domains, 
where U(r) > Uq everywhere (assuming V — a 3 s), at least for the leading term of the asymptotic 
dependence of n s on s (the very use of the continuous model of the random medium suggests that our 
consideration is valid only for s ^> 1). In addition, this estimation could be valid only for clusters with 
Uo 3> a, far away from the percolation threshold, because we consider the distribution of the average field 
in the most compact domain (a sphere). In the vicinity of the percolation threshold clusters typically have 
a fractal-like structure |H]. If this assumption is true, then for the dipolar- like Gaussian field 



U, 



2 



B^s 1 ' 3 , B d = — = 0.34... (5) 



Us * UqsVs ° XP V ~ a a 2 " J' ~ a ~ 4A(36tt)i/3 

We compared (O with the simulation data and found that it provides a good approximation for the true 
cluster numbers (see Fig. [T}. In addition, this approach gives a well-known exact asymptotics (9] for the 
uncorrelated Gaussian field. 
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Fig. 1 Cluster numbers ?i s for the dipolar-like Gaus- 
sian field. Threshold energy Uo ja varies from 2.75 
to 4.0 (with the step 0.25) from the topmost curve 
downwards (inset). In proper coordinates all curves 
approximately collapse to the uniform straight line 
with the slope —0.31 in good agreement with (|5}- 



Fig. 2 Distribution of average squared field Uq in 
the dipolar glass for Ro equals to 5a, 10a, and 20a 
(circles, squares, and crosses, correspondingly). Inset 
shows the plot of the low Uq tails. 



Distribution of the average field in the domain for the Gaussian random field has a Gaussian form, 
irrespectively to the spatial correlation of the Gaussian field. This is not so for the distribution of the 
average squared field 

Qv(Etf) = ( S [yJ drV 2 (r)f v (?) - U^j \ . (6) 
For the uncorrected Gaussian field this distribution, indeed, has a Gaussian form for s ^> 1 

Qv{u ^ ■ = (zi) 1/; 2 ° xp 

while for the correlated Gaussian field this is not so. For example, for the dipolar glass the distribution is 
highly asymmetric (see Fig. |2). For Uq > a 2 the tail of the distribution approximately follows the law 
In Qv oc — (Uq — a 2 ), while for Uq < a 2 the corresponding law is In Qv oc (Uq — a 2 ) 3 . This difference 
directly follows from the strong spatial correlation of U(r): inner parts of clusters provide contribution to 
Qv for Uq > a 2 , while the boundaries between clusters mostly contribute to the tail Uq < a 2 . 



s(ug-*'y 



4cr 4 



(7) 



3 Energetic disorder at the interface: how it differs from the bulk 

Electrostatic nature of the energetic disorder in organic materials has a very nontrivial consequence. Statis- 
tical properties of such disorder at the interface with other material may change dramatically without any 
variation of the spatial or chemical structure of the organic material. 

Let us consider the interface of organic glass with the metal electrode. The electrostatic energetic disor- 
der is directly proportional to the disorder in the spatial distribution of electrostatic potential U (f) = eip(r), 
generated by randomly situated and oriented dipoles or quadrupoles. In organic layers sandwiched between 
conducting electrodes distribution of (p(r) must obey a boundary condition at the electrode surface: at this 
surface the potential should be a constant. Thus, at the electrode surface there is no energetic disorder 
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at all, irrespectively to how disordered is the material in the bulk. This means that the magnitude of the 
dipolar or quadrupolar disorder increases while going away from the interface, asymptotically reaching its 
bulk value o\>. Magnitude of the dipolar disorder at the interface is decreased as (here z is the distance 
from the interface plane) iflOl 



°\z) « o\ [l - g (l - e- 2 */«°) 



a = Aa. (8) 



Asymptotics of the correlation function C(zi, z 2 ,pi — p 2 ) = (U(zipi)U (z 2 , p 2 )) for |ri — r 2 | S> a can 
be calculated using the approach, suggested in ifTTl 

C(z 1 ,z 2 ,p) = a 2 b a l 1 -— = J - ), ( 9 ) 

\ + (ii - z 2 y ^p2 + ( Zl + z 2 ) 2 J 

Here p is a 2D vector oriented along the interface plane. Hence, at the interface the dipolar glass is much 
less correlated in comparison with the bulk: C(zi, z 2 , p) oc ziz 2 /p 3 for p 3> zi,z 2 , while in the bulk 
C(r) oc 1/r, and clusters are elongated perpendicular to the interface. 

Parameter K for a sphere, located at the distance z from the interface, can be estimated by a direct 
generalization of (01 to the case of spatially non-uniform disorder and has a form -Kbuik — SK with 



SK = l67rXao.Ro 



J™ cxp (-2^) 9 2 (k), g(k) = J dtVT^Ji (kVT^) cosh(fci), 



2 v 2 

5K^a b a — , R <C z. 

Here Ji (x) is a Bessel function. Magnitude of K at the interface is smaller than in the bulk, reflecting the 
decrease of the spatial correlation. 



4 Conclusion 

In this paper we discussed the cluster properties of the electrostatic energetic disorder in organic glasses 
with the emphasis on the properties of polar materials built by molecules having permanent dipole mo- 
ments. We found that for this disorder the distribution of the average energy could serve as a good esti- 
mation for the true cluster numbers per lattice site for large "deep" clusters, where s 3> 1 and the thresh- 
old energy Uq is significantly greater than the rms disorder a. Comparison of the analytical results for 
Pv{Uo), calculated for a spherical domain, and computer simulation data for n s supports the hypothesis 
that Pv(Uo) provides the exact leading asymptotic term for n s . We have also shown that the electrostatic 
energetic disorder in organic semiconductors changes dramatically in the neighborhood of a metal elec- 
trode. This phenomenon is of principal importance for the description of the charge injection in organic 
devices. 
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